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Abstract. The current discrepancy of theory and experiment observed recently in muonic hydrogen ne- 
cessitates a reinvestigation of all corrections to contribute to the Lamb shift in muonic hydrogen (/iH), 
muonic deuterium (^D), the muonic '^He ion (denoted here as /^^He"*"), as well as in the muonic ^He ion 
(^t''He+). Here, we choose a semi-analytic approach and evaluate a number of higher-order corrections 
to vacuum polarization (VP) semi-analytically, while remaining integrals over the spectral density of VP 
are performed numerically. We obtain semi-analytic results for the second-order correction, and for the 
relativistic correction to VP. The self-energy correction to VP is calculated, including the perturbations of 
the Bethe logarithms by vacuum polarization. Subleading logarithmic terms in the radiative-recoil correc- 
tion to the 2S-2P Lamb shift of order ci(Za)^ \n{Za) / {nif^ mjv) are also obtained. All calculations are 
nonperturbative in the mass ratio of orbiting particle and nucleus. 
PACS: 12.20.Ds, 36.10.Ee, 14.20.Dh, 31.30.jf, 31.30.jr 



1 Introduction 

In simple muonic bound systems such as muonic hy- 
drogen and deuterium (/iH and //D) and muonic helium 
ions {ji^He^ and /i''He+), the mass ratio of the orbiting 
particle to the mass of the nucleus is larger than the fine- 
structure constant a. The dominant radiative correction in 
these systems is given by electronic vacuum polarization, 
which screens the proton charge on a distance scale of the 
order of the electron Compton wavelength. Here, by "elec- 
tronic" vacuum polarization, we refer to the modification 
' of the photon propagator due to the creation and anni- 
. hilation of electron-positron pairs. Due to the dominance 
' of the vacuum-polarization energy shift (modification of 
I the Coulomb force law at small distances), the 2S level is 
' energetically lower in muonic systems as compared to the 
' 2Pi/2, reversing the ordering found in the hydrogen atom. 
A characteristic property of simple bound muonic sys- 
tems is the mass ratio 



for muonic hydrogen and deuterium, and 
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for muonic helium ions, where the latest recommended 
values of the masses have been used [T] . We here denote 
the masses of the helion nucleus and of the alpha particle 
as TOHe3 and mne* , respectively, in contrast to the muonic 
ions themselves, which we denote as '^He''' and ''^He'*'. 

For heavy muonic ions and atoms, the atomic binding- 
strength parameter is Za, where Z is the nuclear charge 
number and the mass ratio ^at of muon and atomic nucleus 
fulfills the inequality 



Cjv = — < < 1 
rriN 



(3) 



C.N = 

TOat 



of orbiting particle and nucleus, which reads as 



ip = — = 0.112609. 



= — ^ = 0.0563327 . 



(1) 



(2a) 



In heavy muonic ions, the external-field approximation 
(Dirac-Coulomb equation) gives an excellent result. The 
Dirac formalism takes relativistic effects (parameterized 
by Za) into account to all orders. For a system like 
positronium (e'^e^) or true muonium (^+^^, Refs. [2H1]), 
the situation is opposite, 



a = Za « = i = 0(1) 



(4) 
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(2b) each muon is the "nucleus for the other one." In these sys- 
tems, the Breit Hamiltonian is adequate [5]. The (static) 
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Brcit Hamiltonian is exact in the muon-nuclcus mass ratio 
but pcrturbative in Za. For muonic hydrogen, deuterium 
and muonic helium ions, the mass ratio ^at is larger than 
the atomic binding strength parameter Za. These systems 
therefore "lean" more toward the situation encountered in 
positronium and true muonium than toward heavy muonic 
ions. It is thus preferable to treat the reduced-mass depen- 
dence of the corrections exactly (wherever possible). For 
completeness, we would like to stress here that yuH is the 
only muonic atom studied so far (Ref. [6]) where the mass 
ratio is larger than the atomic binding strength parame- 
ter. In all other heavy muonic ions studied primarily in the 
1970s and 1980s (for theoretical overviews see Refs. [ZHH]), 
the binding parameter Za is much larger than the mass 
ratio. 

We proceed as follows. In Sec. [U we evaluate second- 
order as well as relativistic corrections to vacuum po- 
larization using our semi-analytic approach. The latter 
have been the subject of a recent paper [3]. In Sec. |21 
the muon self-energy corrections to vacuum polarization 
are evaluated, taking into account the shift of the Bethe 
logarithms due to vacuum polarization. Finally, radiative- 
recoil corrections are treated in Sec. and subleading 
single logarithmic terms are calculated, supplementing re- 
cent work [10| . Three appendices El— [Q complement the 
paper. We use natural units with = c = eo = 1. 



2 Higher-Order Corrections to VP 
2.1 Second-Order Correction to VP 



where is the electron mass and p is a dimensionless 
spectral parameter for the vacuum polarization. It is also 
useful to define the linear operator A', 
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The reduced mass of the system is written as 
rriN 



P 



ruM 



/3, 



N 



Za fj. 



(6) 



(7) 



The ratio of the muonic Bohr radius to the electron Comp- 
ton wavelength is given by the ratios [llj 



/3p = 0.7373836 . 
= 0.7000861 . 
;3hc3 = 0.3438429 . 
Pn^i = 0.3407691 . 



(8a) 
(8b) 
(8c) 
(8d) 



Using formulas for the reduced Green function from Ap- 
pendix [B] and for analytic integrals from Appendix [CI it 
is possible to analytically evaluate the following matrix el- 
ement, which describes the second-order perturbation due 
to VP. For the 2Pi/2~2S'i/2 Lamb shift difference 



Wvp(Ai;r) 



i'vp(Ai;r) 



1 



(E-H) 



7Wvp(A2;r) 



2P 



2S 



1 



7Wvp(A2;r) 



25 



(9) 



For the calculation of VP effects in muonic systems, it 
is convenient to define the "massive" Coulomb potential ^here the prime denotes the reduced Green function, we 



Wvp(A;r) 



Za 



find 



-Ar 



X = X{p)=meP, (5) 



W2). ^ .7 ^2 / [I + PnIPi + p-lT^ Q 0% [pl + pl + i^N PlP2n ( {I + Pl){l + Pn P2y 

IP^,P2) ^ /^l i2(l + /3^pi)5(l + /3^P2)5 (l + /3ivPir(l + /3ivP2)4 ""l l + /3w(pi+P2) . 



(10) 



Here, Q = Q{Pn', Pi, P2) is a polynomial in the three arguments, symmetric in pi and p2 and rather compact, 
Q = -'ipl + l^Plpt + ^APIpI + %P%pI - 2ApNpip2 - 3/3%pIp2 + 99f3%ptp2 + 111/3% plp2 + 33/3% plp2 
- I^PIpIpI - I^^PIpIpI ~ 33p%p\pl + 39P%p\pI -f I2p%plpl ~ U2p%plpl - &2p%p\pl 
+ M/3%p\pI + 3(!>PIp\pI + lAf3%p\pt + 120/?^pfp^ ^ up^plp^ - l2/3%p\pl + (pi o P2). (U) 



The energy shift is given as where Ki and K2 are the generalizations of the opera- 

tor in Eqs. (jB]) to integration variables pi and p2- Our 
semi-analytic approach allows us to evaluate the remain- 
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Fig. 1. Feynman diagrams for tlie relativistic correction to 
vacuum polarization in the two-body system of a muon and a 
nucleus. The electron-positron pair in the loop is denoted by 
the symbol e. Diagram (a) is the Coulomb photon exchange, 
given by the 00-component of the "massive" photon propa- 
gator, whereas diagram (b) is the magnetic exchange, corre- 
sponding to the spatial components of the "massive" photon 
propagator. 
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(b) 



Fig. 2. Feynman diagrams for the self-energy correction to 
vacuum polarization in a simple muonic bound system. Dia- 
gram (a) represents the vacuum-polarization insertion into the 
exchanged photon in the vertex correction, whereas diagram 
(b) represents the wave function correction to the self energy 
due to vacuum polarization. 

ing two-dimensional numerical integrals to essentially ar- 
bitrary accuracy. 

For the muonic systems of interest, the second-order 
shift is found as 





^ 0.150897 mcV, 


(13a) 




= 0.172023 meV, 


(13b) 


Z\£;(2)(/i3He+) 


= 1.677290 meV, 


(13c) 




= 1.707588 mcV. 


(13d) 


For muonic hydrogen, we confirm the entry 
of Ref. [TT], and for /i^Hc+, we confirm the r 


in Eq. (28) 
esults given 



in Eq. (38) and (39) of Ref. [T2]. Our semi-analytic ap- 
proach ehminates any conceivable numerical inaccuracy 
as the reason for (part of) the experimental-theoretical 
disagreement [5]. 

2.2 Relativistic Correction to VP 

The relativistic correction Sw = J2i=i ^"^i to the Breit 
interaction, relevant for the 2ij/2-2S'i/2 Lamb shift, is 
given as the sum of four terms [TOlITT] 



5wo 



SW3 



Za 



_5i_ 



-\r 



(14) 



Xr 



Za e" 



AmuTriN 



r 



l + \r 



.'■rJ 



Swa ~ Za 



Ami 



2m^m.]s[ 



(1 + Ar) 



cr ■ L , 



where (5/ is one for half-integer nuclear spin and zero for 
integer nuclear spin. For the relevant Feynman diagrams, 
see Fig. [TJ The massive Breit Hamiltonian (|2.2p is a gen- 
eralization of the ordinary Breit Hamiltonian briefly re- 
visited in Appendix |21 The relativistic correction to VP 
is the sum of first-order and second-order perturbations, 





- SE^^'^ + 






(15a) 




= {nlj \5w\ n£j) , 




(15b) 




= 2 l^nej 






(15c) 




= K[6M^^ 




(15d) 


5E,p 


= SE^'^ + 


SE^'^ . 




(15e) 



Here, the reference state has principal quantum number n, 
orbital angular momentum £, and total angular momen- 
tum j. For the second-order matrix elements SM'^^^ one 
needs the explicit formulas for the reduced Green function 
given in Appendix [B] The radial integrals can be per- 
formed analytically using formulas given in Appendix [C] 
Finally, the matrix elements JM*^^) and JM^^^ can be eval- 
uated analytically. 

We present explicit results for the matrix elements 
(5M(^) and SM^'^^ [see Eqs. (jl5bp and (fTSc]) ]. assuming 
(5/ = 1, for the muonic bound systems under investigation, 
in terms of the parameters Za, [defined in Eq. ([Ij], and 
(3n [defined in Eq. ©], 
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2 96 (1 + 4Ar)^ (1 + ^jvp)^ 

, {2 [l + il- (1 + /3jv pf (3 + 4/3^p)] + ^nP (1 + g^) [8 + I3n pjU + 8/?jvp)]} 



5M(2)(2Pi) 



:(i + e7v)Mi + /3A^p)' 

(Za)4 ^ - 8{pNP? (2 + + - 8 (1 + /?N p) (3 + 6^n + 2^^) ln(l + p)] 



, (16c) 



96 (1 + gAr)2 (1 + /3jv p)5 

f{f3N,P,^N)^ -13-23gjv-7e^-/3^p (41 + 67ejv + 19C^) , (16d) 
[5(/?jv, P, Civ) - 8 (Pn pY {^n - 1)^ - 8 (1 + ^n? (1 + Pn p) (l + 2(/3^ p)^) ln(l + Pm p)] 



32 (1 + ew)2 (1 + /3^, pY 



(16e) 



g(/3jv, P, Cw) = - 7 - S^AT - 7^^ - 9/3Arp (3 + Cat + S^at) - 2(/3Arp)2(15 - IB^at + 15^^) - 10(/3wp)3(3 - + 3^^). 

(16f) 



This reduces the calculation of 5E^^'^ and SE'^'^'^ to sim- 
ple numerical integrals over the spectral representation of 
the vacuum polarization {K operator). We have not at- 
tempted here to carry out the remaining integrals over p 
analytically because the integral representation allows us 
to evaluate the relativistic correction to VP to essentially 
arbitrary accuracy. 

For the 2Pi/2~2S'i/2 Lamb shift, we indicate the dif- 
ference as Z\£'vp = J£'vp(2Pi/2) - (5£'vp(25i/2), 



Zi£;vp(pH) 


= 0.018759 meV, 


(17a) 


^-Bvp(pD) 


= 0.021781 meV, 


(17b) 


Z\i;vp(p'He+) 


0.509344 meV, 


(17c) 




= 0.521104 meV. 


(17d) 



The relativistic correction to vacuum polarization is repre- 
sented by the tree- level Feynman diagrams in Fig.[TJ which 
represent vacuum-polarization insertions in the Coulomb 
and magnetic exchange, with relativistic wave functions in 
the in and out states. Our results are in agreement with 
those recently reported in Ref. [TU]. This concludes our 
semi-analytic rederivation of the relativistic corrections to 
vacuum polarization, with proper account of the reduced- 
mass dependence. 



3 Self-Energy Correction to VP 



for the Schrodinger-Coulomb problem, evaluated on ex- 
act operators that take vacuum polarization into account. 
Here, we compare the nonperturbative to a pcrturbative 
treatment of the vacuum-polarization potential Vvp. The 
Uehling potential is added to the Schrodinger Hamilto- 
nian by the replacement V V + Vvp • The effect of high- 
energy virtual photons in the self-energy loops given in 
Fig. [5] can thus be expressed in terms of the Dirac Pi 
form factor acting on the vacuum polarization potential 
Vvp . When rewritten in terms of the noncovariant photon 
energy cutoff e, which is a convenient overlapping param- 
eter in Lamb shift calculations |14j . we have in first order 
in Kp, 



5E 



a 



H 



In 



\2e J 9 



+ 2 ( n£, 



vp 



E-H 




(18) 



This is equivalent to expanding Eq. (3.7) of Ref. [T^ to 
first order in the vacuum-polarization potential. The cor- 
rection 6Eh to the high-energy part is expressed in terms 
of a parameter Vgi , 



In Sec. 3.2 of Ref. [13], the self-energy correction to 
vacuum polarization has been discussed. In the cited ref- 
erence, a partially nonperturbative approach has been 
chosen in order to evaluate certain matrix elements in a 
nonperturbative framework, using exact wave functions 



SE 



H 



^6i^n(^)+^^. (19) 



10' 
9" 
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Wc find, using techniques similar to those employed in 
Seel 



logarithm due to the Uehling potential. Numerically, we 
find for 2Fi/2, 



F6i(2Pi/2;AiH) = 


- 0.02327, 


(20a) 




- 0.014559 , 


(25a) 


V6i(2Pi/2;AiD) = 


- 0.02453 , 


(20b) 


Lfin(2Pi /9: uD) = 


- 0.015446 , 


(25b) 


V'6i(2Pi/2;M^He+) = 


- 0.04345 , 


(20c) 


Lfin(2Pi/9: w^He+) = 


- 0.032293 , 


(25c) 


*^6l(,^A/2i M j — 


— U.U4oDy . 


\i\j(X ) 


7" f C\ T~t 4-TT -!- \ 

-^6o(2Pi/2;Ai He+) = 


— 0.032573 , 


/or 1 \ 

(25d) 


For the 2iS'i/2 state, the results are 






whereas for 2Si/2, 






^61 (2^1/2 ;mH) = 


3.08601, 


(21a) 


Leo (2^1/2 ;AiH) 


= 11.176, 


(26a) 


1/61 (2^1/2 ;a^d) = 


3.18785, 


(21b) 


L6o(2S'i/2;/iD) 


11.464, 


(26b) 


^61 (2^1/2; M3He+) = 


4.71872, 


(21c) 


L6o(25i/2;/i'Hc+) 


= 15.640, 


(26c) 


^61 (2^1/2; A^'He+) = 


4.73968 . 


(21d) 


L6o(25i/2;M'He+) 


= 15.696. 


(26d) 



The results for /iH confirm the entries in Eq. (3.8) 
of Ref. |13| . where a nonperturbative approach in the 
vacuum-polarization potential was employed. The correc- 
tion due to the anomalous magnetic moment of the elec- 
tron reads as 



The total self-energy vacuum-polarization correction to 
the 2P-2S' Lamb shift then is 



a\Zaf ^3 



SE.^p = SEh + SEm 



SEr = 



(27) 



SEm = 



IdV, 
r dr 



^(<T.L) 



In 



10 



AVei + AMe, 



4 

^ALeo 



+ 2 ( ni-i 



1 



where 



ne,, 



(28a) 



60 • 



(22) 



The matrix element Meo is nonvanishing for P states. 



M6o(2Pi/2;A'H) = 
M6o(2Pi/2;mD) = 
M6o(2Pi/2;M'He+) = 
Af6o(2Pi/2;M'He+) = 



0.04276 , 
0.04656, 
0.13396, 
0.13556. 



(23a) 
(23b) 
(23c) 
(23d) 



ZiHi =yei(2Pi/2)- -^^61 (2^1/2), 
Z\Meo = il/6o(2Pi/2) - il/6o(2^i/2) , (28b) 

Z\Leo = L6o(2Pi/2) - L6o(25i/2) . (28c) 

The final numerical values for the contributions to the 
2Pi/2^2S'i/2 Lamb shift are given as 



For S states, we have M6o(25i/2) = in view of angular 
symmetry. For the 2P state, we here take the opportunity 
to correct the result given in Eq. (3.10) of Ref. [13]. The 
correction has negligible infiuence on the numerical result 
reported below in Eq. (|29p for muonic hydrogen. 
The low-energy part is conveniently expressed as 



Z\S,vp(/iH) = 


- 0.00254 meV, 


(29a) 


Z\P,vp(mD) - 


- 0.00306 mcV, 


(29b) 


Z\P,vp(M'He+) = 


- 0.06269 mcV, 


(29c) 


Z\P,vp(M^He+) = 


- 0.06462 mcV. 


(29d) 



8Et 



4 ,,3 r 



1^61 In 



{Za)^H 



■Lfi 



(24) 



where the Vgi term leads to a cancellation of the e parame- 
ter. The coefficient Leo gives the modification of the Bethe 



4 Recoil Correction to VP 

4.1 Formulation of Radiative Recoil 

The first genuine two-body energy correction beyond 
the frequency-independent part of the Breit interaction in- 
volves a photon-frequency dependent transverse exchange 
and a relativistic two-photon exchange [T5]. For the IS' 
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(a) 




(c) 

Fig. 3. Typical Feynman diagrams for the recoil correction. 
Diagram (a) is the frequency-dependent part of the Breit in- 
teraction, (b) is the seagull exchange graph, and (c) and (d) 
are two-photon exchange graphs. 



state, the corresponding energy shift was calculated in 
Ref. |16| . which is why the recoil correction is commonly 
referred to as the Salpeter correction, and the result was 
generalized to an arbitrary excited state in Ref. [T7] . The 
correction reads, for an individual state with principal 
quantum number n and angular momentum £, 



2 . 

3 5. 



In 



1 \ 8 (5,0 

— — in fcn 

Za i 9 



_7a^ 
3 

where 



25, 



m 



N 



2 1 (""^N 

m,, in 



'TV 



In 



In 



+ 1 + _L + ^(^ri) + JE 
2n 



(30) 



(31) 



and 'F{n) = X]fe=i ^"^ denotes the logarithmic derivative 
of the Gamma function. 

For reference, we evaluate the recoil correction for the 
2Pi/2~2Si/2 Lamb shift, for the muonic systems under 
investigation, 





- 0.044971 meV, 


(32a) 


AEnifiJ^) = 


- 0.026561 meV, 


(32b) 


AEii{p^Rc+) = 


- 0.558107 meV, 


(32c) 


AEr{^,^Rc+) ^ 


- 0.433032 meV. 


(32d) 



We thus confirm the result indicated in Eq. (66) of 
Ref. [12] and the results given in Tables 3-6 of Ref. [H] . 

According to the derivation presented in Ref. [TS] , the 
recoil correction is formulated in terms of four corrections. 




^0^ 



p p 

(b) 





p 


p 


p 




P 


(c) 






(d) 

















Fig. 4. The six Feynman diagrams for the vacuum- 
polarization correction to recoil (radiative-recoil correction) 
of relative order a, in muonic and antiprotonic systems. The 
electron-positron pair in the loop is denoted by the symbol e. 
Coulomb photons are denoted by dashed lines. 



Two of these are associated with the photon-frequency 
dependent part of the Breit interaction [see Fig. ^a)]. 
The other two are associated with two-photon exchange 
[sec Fig. [3l[b), (c) and (d)]. The frequency-dependent one- 
photon exchange in Fig. ^a.) leads to a low-energy part 
El where the photon frequency fulfills uj <^ Za fi and 
therefore is small against the orbiting particle momen- 
tum, and a middle-energy part Em, where lu ^ Zaji. 
The two parts arc separated by an overlapping parame- 
ter. The frequency-dependent part of the Breit interaction 
is described by El and Em- 

For low photon frequencies, the two-photon exchange 
is described by the seagull graph in Fig. [31[b) . The seagull 
exchange energy correction Es effectively represents the 
"low-energy part" of the two-photon exchange, whereas 
the high-energy part Eh describes a local operator, with 
both virtual photon momenta and the momenta of the 
constituent particles being large [see Figs.|3Uc) and (d)]. It 
contributes only for S states. A second overlapping param- 
eter cancels between Es and Eh- The seagull exchange di- 
agram in Fig.|3Jb) exists in nonrelativistic QED (NRQED) 
as opposed to fully relativistic QED, where two-photon 
emission out of the same vertex is forbidden. Indeed, the 
seagull diagram follows from the two-photon exchange di- 
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agram of rclativistic QED (two-photon emission) in the 
hmit of soft photons and negative-energy virtual fermion 
states (the fermion lines would have a "Z" shape in time- 
ordered perturbation theory) . 

The radiative-recoil correction of order 
a{Za)^fj,^/{m^mN) is the perturbation of the recoil 
correction of order [Za)^ jj,^ / {m^mN) by an ordcr-a 
vacuum polarization, and is calculated here for simple 
muonic bound systems. Because the correction is nu- 
merically small, we restrict our attention to the leading 
logarithms. The contributing diagrams at relative order a 
are given in Fig. |4l One might expect a seventh diagram 
to contribute, with a VP correction to a Coulomb photon 
exchange in the seagull graph, where the seagull vertices 
connect to the muon and proton/deuteron lines on op- 
posite sides of the Coulomb/Uehling exchange. However, 
the seventh diagram only contributes at relative order 
a{Za) and is not considered here. 



4.2 Logarithm from Middle- and Low-Energy Part 



As evident from the derivation in Ref. [T5], the loga- 
rithmic term in Eq. (jSOj can be broken down as follows. 



E 



{Zaffi'^ (8 



log _ 

^ ■nm^misin'^ 



- 2 1 (5,0 In 



1 

Za 



(33) 



where the term with the prefactor "8/3" comes from the 
sum of the middle- and low-energy parts. The first loga- 
rithmic term in the recoil correction (due to the middle- 
and low-energy parts) can thus be expressed as a matrix 
element of the Laplacian of the Coulomb potential V ^ 



For the 2Pi/2-2S'i/2 Lamb shift, indicated by the prefix 
Z\, the correction evaluates to 





- 0.000505 meV, 


(36a) 


AEL{^i^) = 


- 0.000295 meV, 


(36b) 


AEL{fi^lic+) = 


- 0.005724 meV, 


(36c) 


AEL{fi^llc+) = 


- 0.004431 meV. 


(36d) 



This concludes the treatment of the logarithmic terms 
in the radiative-recoil correction due to the diagrams in 
Fig. ^a.) and (b). It is not the leading logarithmic term 
in the radiative-recoil correction. 



4.3 Logarithm from Seagull and High-Energy Part 

The mechanism for the generation of logarithmic terms 
in the seagull part is different from the middle- and low- 
energy parts and cannot be traced to a perturbative mod- 
ification of the matrix element in the term ([M]). For the 
first logarithmic contribution, it is sufficient to consider 
the seagull exchange diagram given in Fig.|4l^c); the high- 
energy terms from Fig. HKe) and (f) cancel an intermedi- 
ate overlapping parameter and their contribution is not 
needed in the evaluation of the logarithmic terms. The 
seagull term, with a vacuum polarization insertion in the 
exchanged photon, is given as 



SEs 



K 



(27r)3 (27r)3 ui fcz wi 



(37) 



hi U 



i (fci+fc2)-'r 



E, = . {ni, I VV| nti) In f ^ ) . (34) 



37r TTiu totv 



Perturbing the matrix element by the vacuum- 
polarization potential, we obtain the following logarithmic 
(L) term in the radiative-recoil correction. 




Ho 



vp 




a 2{Zafj^ ( 1 

TT TT mj, TL'^ \ Za 



(35) 



This result is spin-independent and in full agreement with 
Eq. (5.165) of a recent unpublished work (Ref. [IH]), where 
for P states, nonlogarithmic terms are calculated as well 
(we only consider the leading logarithms here) . Results for 
Vfji can be found in Eq. (j20ap . 



where oji — \/k\ + \^ = ^Jk\ + {irie p)^ is the photon 
frequency for massive photons. From this integral, using 
a cutoff A as in Eq. (27) of Ref. [T3] as an overlapping 
parameter that separates the integration region from the 
high-energy part, we extract the following logarithmic cor- 
rection terms, which contribute only for 5* states. 



4Q(Za)V3 
oEs = T-^ ^ Ou 

X [- In^ (4Za/3 



(38) 



ln(2/3Ar) \n{AZaPl)] . 



The double logarithmic term is formally leading (we here 
confirm the result of Ref. [IH])- This correction evaluates 
to 



AEs{^m) 


= 0.000414 meV , 


(39a) 


AEsifiD) 


= 0.000251 meV, 


(39b) 


AEs{fi^}ie+) 


= 0.006648 meV, 


(39c) 


AEs{fi^}ic+) 


0.005174 meV. 


(39d) 
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There is a second correction for S states. It is given by the 
wave function correction to the leading seagull logarithm 



Es 



2{Zaf 



(40) 



where 0„f is the nonrelativistic Schrodinger wave func- 
tion. We note that this term cannot be traced to the ex- 
pectation value of the Laplacian of the Coulomb potential, 
(V^T^), and therefore the formalism leading to the result 
in Eq. psp is not applicable. The wave function correction 
due to the potential Wvp is 



E -H 



•nSI 



(41) 



At the origin, the correction can be expressed as (25" state) 
(ZaM)5/2 



4V2^(A + Za^if 
+A{Za^Lf\^ + ll{Za^lf\ + 3(aAi)^ 
^ {Za^,)V^{2\^ + {a^,f) 



27r(A + Za^iY 



Za jji 



(42) 



The corresponding logarithmic energy correction for an 
nS state due to the wave function correction reads as 



5E 



w = -^^^i^[20„,(O),5<^„,(O)]lnf^) (43) 
m^^raj^ \Za, 



for an individual state. A numerical evaluation for the 
Lamb shift (2Pi/2~2S'i/2 energy difference) leads to the 
results (again denoted by the prefix A instead of 5), 

AEw (/iH) = 0.000228 meV , (44a) 

AEw (AiD) = 0.000138 meV , (44b) 

AEw{^i^^c+) = 0.004017 meV, (44c) 

AEw{^i'^^c+) = 0.003123 meV, (44d) 

which is of the expected magnitude. The correction AEw 
corresponds to the diagram in Fig. IDJd). The seagull 
and high-energy parts do not generate any logarithmic 
radiative- recoil corrections for P states. 



4.4 Total Logarithmic Radiative-Recoil Correction 

For the Lamb shift, the total logarithmic radiative- 
recoil correction is given as the sum 



AE 



RR 



AEt. + AEs + AE, 



w ■ 



(45) 



It evaluates to 

AERR{^i}^) = 0.000136 meV, (46a) 

AEiiRinB) = 0.000093 meV , (46b) 

AEiiii{ti^}ic+) = 0.004941 meV, (46c) 

AEiiii{ti^}ic+) = 0.003867meV. (46d) 

The results are numerically small and suppressed with re- 
spect to the leading recoil correction given in Eq. p2[) by 
a factor a. 

There is some numerical cancellation between the lead- 
ing squared logarithm given by the first term on the right- 
hand side of Eq. ((38)l and the subleading single logarithms. 

5 Conclusions 

The conclusions of this paper are twofold. The first 
conclusion is that a number of nontrivial higher-order cor- 
rections to the Lamb shift in simple muonic systems (/iH, 
fiD, ^"^He^, //'*IIe+) can be evaluated semi-analytically. 
Certain well-defined, remaining one- or two-dimensional 
numerical integrals can easily be evaluated to essentially 
arbitrary accuracy. The second-order VP shift (Sec. 12. ip 
has been investigated numerically in a number of previ- 
ous works |11H13| , and we here confirm these results using 
our semi-analytic approach [see Eq. and complement 
the literature with a result for ^'^IIe+. The second-order 
correction is otherwise easy to evaluate numerically; the 
only advantage of our approach is that the integral rep- 
resentation easily allows us to evaluate the correction to 
essentially arbitrary accuracy, yielding an additional con- 
firmation for the literature values. The relativistie cor- 
rection to VP (with a proper account of the reduced- 
mass dependence) is analyzed in Sec. 12.21 The calcula- 
tion is reduced to one-dimensional parametric integrals 
[see Eqs. ([TS]) and P^ ]. These expressions yield additional 
evidence for the results given in Eq. p7|l . We have not 
attempted to find analytic expressions for the remaining 
one- and two-dimensional integrations over the spectral 
parameter p of the vacuum polarization, but note that 
such calculations may be possible j20,.21]; we leave this 
problem for future investigations. 

The second point of the paper is the calculation of two 
nontrivial higher-order corrections to vacuum polarization 
in simple muonic bound systems: namely, the self-energy 
and recoil corrections to VP (see Sees. [3] and |4]). For the 
self-energy correction to VP, the final numerical results 
are given in Eq. (|29p . Numerically, we confirm results re- 
ported previously in Ref. [13] for fiH. The same analytic 
techniques are used as those employed in Sec. [2] but in- 
termediate results are suppressed. In Eq. (j23ap . we take 
the opportunity to correct |13| a computational error in 
the evaluation of the self-energy correction for the matrix 
element A/go which describes the vacuum-polarization cor- 
rection to the electron anomalous magnetic moment term 
in the Lamb shift (the correction has negligible effect on 
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the numerical value of the final result). Our results include 
the correction to the Bethe logarithm due to vacuum po- 
larization [see Eq. (psj) and (HH)]. 

The radiative-recoil correction is found to be numeri- 
cally small for simple muonic bound systems. Leading and 
subleading logarithmic terms are calculated here. The first 
logarithmic term in the radiative-recoil correction is due 
to the vacuum-polarization correction to the frequency- 
dependent Breit exchange and is given in Eq. ([55]) [see 
Figs.HI^a) and (b)]. From the seagull part, the leading log- 
arithm is given (for S states) by the vacuum-polarization 
insertion into the exchange photon lines [see Figs. Hl^c), 
(e), (f)]. The corresponding double- logarithmic term SEs 
can be expressed analytically, scales as where n is the 
principal quantum number, and can be found in Eq. psp . 
The appearance of the parameter /3n in the result implies 
that the logarithmic coefficient depends on the mass ra- 
tio of muon and electron, as it should. Finally, the wave 
function correction to S states which is nonvanishing, gen- 
erates a third logarithm, SEw, given in Eq. P5|) [see also 
Fig-EUd)]. The total results for the radiative-recoil correc- 
tion (logarithmic terms) for the systems under investiga- 
tion can be found in Eq. (j46p . 

Self-energy corrections to vacuum polarization and 
radiative-recoil corrections are found to be of minor sig- 
nificance, in accordance with simple order-of-magnitude 
estimates. Still, the calculations are necessary in order to 
exclude a conceivable large logarithmic term in the higher- 
order effects as an explanation for (part of) the observed 
spectroscopic discrepancy Also, the calculation of the 
subleading logarithmic terms in the radiative-recoil cor- 
rection clarifies the size of one of the traditionally most 
elusive corrections for two-body bound systems, in the 
case of a bound muon. Experiments on these systems are 
ongoing at PSI (Paul-Scherrer-Institute, Villigen). 



shift, the relevant terms in the Breit Hamiltonian read 



5H = 
SH2 = 

5Hi = 



SHi = - 



P 

8m3 



P 



8m% ' 



1 



Za 



2771^ m^v 
Za f 1 



r 



cr ■ L . 



(47) 



where the summation convention is used. Here, Sj = 1 for 
half- integer and Si — for integer nuclear spin (see [12). 
The expectation values of of the Breit Hamiltonian, eval- 
uated for the 2Pi/2^25'i/2 difference read as 



L(2Pi/2-25i/2) 



(Za)V 
48 m% 

{Zafix^ 



(48) 



Si = 0., 



where we recall that 2Pi/2 is energetically higher. The 
shift evaluates to 



L(^H) = 0.05747 meV, 

L(/iD) = 0.06722 meV, 

L(M^He+) = 0.12654 meV, 

L(M^He+) = 0.29518 meV. 
These results confirm entries in Ref. [TOl. 



(49a) 
(49b) 
(49c) 
(49d) 
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A Breit Hamiltonian 



The Breit Hamiltonian describes the fine- and hyper- 
fine splitting in two-body systems in the order (Za)'^, ex- 
act to all orders in the mass ratio. For the 2P-2S Lamb 



B Reduced Green Functions 

We use the reduced Green function G' of the 25" state 
for general radial arguments, 



1 



r2 



, Hq — E2S , 

where R2s{ri,r2) has the following representation, 
afj."^ exp (-i (i?< + P>)) 



(50) 



R2siri,r2) 



4:R< P5 



X [4e^< i?< (P> - 2)P> +Q(i?<,i?>) 

-4i?<(i?< -2)i?>(i?,> -2)P(P<,P>)] . (51) 

Here, i?< = afir^ and i?> — afir^, with r< ~ min(ri,r2), 
and r> = max(ri, r2). The result (j5ip constitutes a confir- 
mation of Eq. (23) of Ref. [TT]. The function F is expressed 
as 

P(P<,P>) Ei(P<) - 27B - ln(i?< i?>). (52a) 
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The exponential integral can be written as 

/■^ e' - 1 
Ei(x) = / dt——+-/E+H\x\) 



= -(P.V.)/ dt— --Ei(-x), (52b) 

for real x (positive or negative), where = 0.577216 . . . 
is the Euler-Mascheroni constant. We denote the principal 
value by the symbol "P.V." . The polynomial Q in Eq. (|51l) 
reads 

g(i?<, i?>) = -8 (i?< + i?>) + 26i?< i?> (i?< + i?>) 
+ R% Rl (i?< + i?>) + 4 (i?^ - 3R<R> + RD 
- i?< i?> {2R\ + 23i?<i?> + 2Rl) . (53) 

Using the integrals listed in Appendix \^ it is actually 
possible to calculate the wave function correction 



1 



EnS — He 



Vvp\4>ns) 



(54) 



analytically, where Wvp is given in Eq. ([5]). 

C Integrals 

In order to carry out the radial integrations in the 
calculation of the wave function perturbation 



1 



Eq — Hq 



(55) 



for 25* and 2P, one needs integrals with integration do- 
mains r € (0,s), and r G (s, oo), due to radial ordering. 
We give the results for two useful integrals of the first 
category: 



(56b) 



dre""'' ln(r) 

= - [Ei(-as) - 7£ - ln(a) (e"'^-^ + l)] , 
a 

I dre-"'^Ei(&r) 

1 

a 

In the second category, the following expressions are rele- 
vant: 

/ dr e-° ln(r) ^ - [Ei(-a s) ln(a) - Ei(-a s)] 
Js a 



Ei ((5 - a) s) - e-" " Ei {b s) + In 
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( 


a — b 


)1 



(57a) 

f°° 1 

/ dr e"'' Ei(6 r) ^ - [Ei ((5 - a) s) - e"'' ^ Ei {b s)] , 

J s 

(57b) 

where we note that the results are much more compact 
than for the first category. Finally, integrals over the in- 



terval r S (0,oo) are needed for the evaluation of second- 
order matrix elements [see Eq. ([9])]. One example is 



dre-'''"Ei(&r) ln(r) = — + — In^ [a{a - b)] 
3 a 2 a 

+ - ln(a — 6) 7_E + In f ^ 
a I V 



7s In(fe) + Li2 



a — b 



(58) 
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